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SUMMARY 

The status of research on the transient behavior of finned coil cross-flow heat exchangers 
using single phase fluids is reviewed. Applications with available analytical or numerical 
solutions are discussed. Investigation of water-to-air type cross-flow finned tube heat exchangers 
is examined through the use of simplified governing equations and an up-wind finite difference 
scheme. The degenerate case of zero air-side capacitance rate is compared with available exact 
solution. Generalization of the numerical model is discussed for application to multi-row multi- 
circuit heat exchangers. 


INTRODUCTION 

Significant activity has occurred over the last 25 years in the study of the dynamic 
behavior of heat exchangers. Analytical and numerical solutions of a variety of heat exchanger 
geometries and applications have appeared. The focus here is specifically with respect to cross- 
flow heat exchangers for application in the HVAC (and automotive) areas; where geometries 
are significantly smaller than found in power generation and process industries. 

In-depth study of the dynamic behavior of such heat exchangers commenced in the mid- 
sixties through research aimed at improving dynamic control capabilities of air heating and 
cooling systems. Numerous simplified first-order transfer function models emerged, such as for 
bare single tubes in cross flow (refs. 1,2), single pass finned tubes, (ref. 3), serpentine single 
row (refs. 4,5,6 to 10), and multi-pass cross-counter and cross-parallel flow (refs. 11,12) heat 
exchangers. Large multi-pass and shell and tube heat exchanger dynamics are discussed in refs. 
13 and 14. These references apply to one fluid mixed, liquid to gas cross flow heat exchangers. 
Most of the results presented are in gain and time constant / frequency response format most 
useful for control engineering. The aim of all of the above mentioned work was to derive 
transfer functions relating outlet temperature response of the primary and secondary fluids to 
changes in flow rate or temperatures of the inlet conditions of primary or secondary fluids. 
Each of the above references provide up to 4 transfer functions relating the responses of the 
outlet temperatures to changes in the inlet flow rates and temperatures. Ref. 15 provides for 6 
transfer functions, accounting for primary and secondary fluid temperatures and primary fluid 
flow rate effects on the outlet temperatures. Ref. 16 provides for 8 transfer functions, including 
the effects of secondary fluid inlet flow rate on the outlet temperatures. Models are also 
available which include closed-loop feed back analysis of heat exchanger transients (ref. 17) and 
heating and cooling system simulation (ref. 18). 

Most of the above works account for the presence of external fins, but do so by 
"lumping" the thermal mass of fins with the tubing, thus neglecting thermal diffusion lags due 
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to conduction in the fins. Work by Waede Hansen and Demandt (ref. 19) and Kabelac (ref. 20) 
explores distributed effects of heat diffusion in the fins on the overall response of finned tube 
heat exchangers. Waede Hansen and Demandt concluded that neglecting heat diffusion effects 
in the fins was reasonable in the frequency range up to 10-15 rads/sec. Kabelac achieved good 
agreement with experiment using a model which incorporated fin diffusion effects, but did not 
separate the effects of the fins from the overall model behavior. 

Many of the more rigorous analytic and numerical approaches for transient heat 
exchanger analysis have been developed more recently and have results presented in temperature 
vs. time format (i.e. the focus is on heat transfer behavior). Some of these studies are more 
applicable to a wider class of crossflow heat exchangers (i.e. plate heat exchangers, rod bundles, 
etc.). Rizika (ref. 21) analytically examined a tube in a constant temperature or insulated 
environment. Evans and Smith (ref. 22) analytically investigated single pass cross flow heat 
exchangers with neither fluid mi xed, bu t did not include th e effects of metal or matrix thermal 
capacitance. Myers et al. (refs. 23,24) developed an approximate integral method for single 
row, single pass geometry (solutions developed for both fluids unmixed). Myers et al. (ref. 25) 
numerically investigated the case of one fluid having an infinite capacitance rate. Jang and 
Wang (ref. 26) numerically investigated the case of one fluid mixed. Terasaka et al. (ref. 27) 
used an approximate method of weighted residuals to investigate the case of both fluids unmixed, 
neglecting metal/matrix capacitance. Spiga and Spiga (ref. 28) utilized the Laplace transform 
technique (and numerical inversion) to investigate the case of both fluids unmixed, with 
metal/matrix capacitance incorporated. Gvozdenac (ref. 29), Spiga and Spiga (ref. 28), Romie 
(ref. 31), and Chen and Chen (ref. 32) used similar Laplace transform approaches to derive 
analytical solutions for gas-to-gas, neither fluid mixed geometries. Yamashita et al. (ref. 33) 
used a finite difference technique to numerically investigate heat exchanger transients with both 
fluids unmixed. Chiang et al. (ref. 34) numerically investigated a cross counter flow automotive 
application. Numerical investigation of multi-row cooling coils with condensation occurring on 
the fins has been studied by Reichert et. al (ref. 35). 

While significant effort in the above works have been applied to developing analytical 
solution techniques and presenting the general transient behavior of the broad class of cross flow 
heat exchangers, less emphasis has been applied to quantifying behavior of particular applications 
and circuiting arrangements. Numerical analysis remains the strongest and most flexible 
technique for examining heat exchanger transients. In the same vein as the work of Myers et 
al. (ref. 23,24, analytical) and Reichert et al. (ref. 35, numerical, cooling coils), the present 
study further examines the analysis, mathematical/numerical description, and behavior of heat 
exchangers for air heating in HVAC. 


MATHEMATICAL DEVELOPMENT 

A length of finned tube exposed to air in crossflow is to be modelled. The tube could 
be coiled in serpentine manner, laying in a plane which is perpendicular to the external fluid 
flow, or it could be a single length of tube in cross flow. Situations where multi-row or 
parallel/counter cross flow circuiting is employed will be discussed later. The following 
assumptions will be applied: 
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1. Incompressible flow, 

2. "Plug flow" model for water. 

3. Convective coefficients constant during transients, equivalent to steady values. 

4. Negligible axial conduction in tube and fluids. 

5. Negligible radiative effects. 

6. Fins/tube are at one unique temperature at a given x location ("lumped thermal mass"). 

7. The temperature difference between the tube/fin surface and the air is the logarithmic 
mean temperature difference. 

8. Negligible thermal storage in the air. 

The primary fluid (water) is treated as mixed (no variation in temperature except in the flow 
direction) and the secondary fluid (air) is treated as unmixed (variations in temperature both in 
the flow direction and normal to the flow direction). 

Energy balances on an elemental unit of fin/tube (m), water (w), and air (a) gives the 
governing equations in the primitive variable form: 

dT — 

K.tf * - W- - T .) * pyA'tf. - rj - 0 (1) 

dT dT 
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As will be shown later, this form may be of most usefulness for studying individual effects of 
physical parameters. The case to be studied here will be for the secondary fluid (air) being 
heated. Thus, temperatures are non-dimensionalized according to (T-TJ/(T wr TJ. Substitution 
and rearranging gives: 
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To reduce the number of variables, dimensionless groupings as follows are helpful: 

X a x/L , j- = tVJL 
K 3 KPUrhj:^ 
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Insertion into the governing equations produces 
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Substitution of Eqn.(9) into Eqns. (7) and (8) yields 
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where f(NJ = 1 -exp(-NJ . Initially, the primary fluid, metal, and secondary fluid temperatures 
are at zero. The primary fluid (water) temperature is subject to a sudden step change in 
temperature. The boundary and initial conditions are then 


0Jx,0) = 9Jx,0) = 0 


( 12 ) 


K(l ,r) - 1 (13) 

The dimensionless groupings are the minimum number required to describe the problem. While 
N W =NTU W and N^NTU, have the same form as the conventional "Number of transfer units", 


they are not actually NTU values, 
conventional NTU (ref. 36) can be 

Conversion of these dimensionless 
accomplished as follows: 
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N, and N w are chosen instead of NTU because they appear naturally in the governing equations, 
and effectively separate primary and secondary side fluid effects. N t dictates how close the air 
temperature comes to the metal temperature (see eqn.6)„ N w contains the water flow length L 
in it, while N, does not explicitly contain the water flow length. This is a result of the energy 
equation for the air being lumped and not distributed; thus no characteristic length appears 
explicitly on the air side. Use of the steady state NTU is not straight forward in parametric 
studies where the definition can vary. N, and N w are also used by Spiga and Spiga (ref. 28) and 
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Chen and Chen (ref. 32). As will be shown, even the non-dimensional parameters selected do 
not allow for a perfect separation of all fluid and heat transfer effects. For example, variations 
in air side resistance affect both the N, as well as H,. This duplicity cannot be avoided however, 
even with a different selection of non-dimensional parameters. 


NUMERICAL TECHNIQUE 

A forward time backward space finite difference scheme was chosen for the numerical 
calculations. The explicit formulation was chosen as a "first cut" approach. Eqns. (10) and (11) 
become 
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A stability criteria arising from Eqn.(17) is seen to be 
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N w in Eqn.(18) arises from the convection term in Eqn.(ll). In a purely advective problem, 
N w =0 and Ar=Ax would result (a desirable result, since in dimensional form it is equivalent 
to At-V w =Ax, which implies that a fluid front properly arrives at a spatial grid point precisely 
at the time dictated by the fluid velocity.) In the present work the numerical constraint on At 
is Ar<Ax, which implies a fluid front will appear at a given grid point prior to what is 
physically dictated by the fluid velocity. Note that by writing Ar=</>-Ax, where 0<4>< 1, for 
best representation of the physics, we want <f> * 1. We can then methodically select grid spacing 
according to 


Ax - (19) 

<t>N w 

and as <t> approaches 1, At approaches Ax- Computational effort increases significantly as this 
is done. 


RESULTS AND DISCUSSION 

A parametric study is presented using the above parameter definitions over ranges of the 
parameters found to cover the design and operation of several manufactured air heating coils as 
determined by a survey of manufacturer’s catalogs and designs reported previously in the 
literature. 


99 



An exact solution for the problem at hand is not available. However, Myers et al. (ref. 37) 
provide an analytic solution for the degenerate case of zero air side capacitance rate. This 
idealized situation was simulated in the present case by forcing N, to be very large. 

Fig. 1 provides a comparison between the numerical model and the exact solution for the 
case of negligible air side capacitance. Fig. 1 illustrates a severe test of the numerical model, 
that of correctly predicting the step response of the primary fluid (water) at the outlet, at a water 
flush time of one. The comparison is seen to improve as the value of <f> in Eqn. 19 approaches 
1.0 (i.e. decreasing grid and time spacings). A reasonable trade-off between accuracy required 
for this study and computational effort was to keep <t> in the neighborhood of 0.995. Figure 2 
illustrates that the 4> criteria has a less dramatic effect on the secondary side fluid (air). For a 
fixed <t> less than 1.0, error in temperature is approximately constant versus time, although the 
magnitude of temperature increases with time and the error appears larger at longer times. The 
ability to control numerical error is demonstrated; however, since numerically the 4> ^ 1 implies 
some error, the explicit formulation is not particularly attractive on this point and future studies 
should explore implicit techniques. 

Figures 3 and 4 reflect the effects of the metal to water capacitance ratio on outlet water and 
air temperatures, respectively. The value of C r does not effect the steady state values of the 
dimensionless temperatures, but has a strong influence on the transient behavior as expected. 
Response times can be compared directly between each of the curves on either figure, with 
greater response times obtained for larger C r (larger wall capacitance). 

Transient heat exchanger efficiency is a useful parameter which compares the energy leaving 
the heat exchanger via the air stream to the energy that is given up by the water stream, at a 
given instant during the transient. At steady state, the efficiency will equal one. For heating 
of the air, we define the transient efficiency as: 


(mc\ (T - T \ 

Transient Efficiency , e - — 111 \ - ± 

KL(n, - rj) 

Figure 5 reflects the transient efficiency versus time for the case of varying capacitance ratio. 
Discontinuities in the curves appear at a dimensionless time of one due to the fact that the water 
outlet temperature changes abruptly at one flush time. While discontinuities are not an attractive 
feature for performance diagrams, this is a reflection of what is really occurring. The transient 
efficiency plot also allows us to compare overall heat exchanger response times for different 
values of the parameter being studied. 

Figures 6 and 7 show the dependence of outlet water and air temperature on N„ respectively. 
One can view the N, variation as a variation in the air flow rate, with smaller N, values 
reflecting larger air mass flow rates (with the air side resistance artificially held fixed). As N, 
decreases, the water outlet temperature, and hence the air outlet temperature, are seen to 
decrease. Figure 8 depicts the transient efficiency from which overall response times may be 
compared. Comparisons of Figs. 6,7, and 8 demonstrate the utility of presenting both outlet 
temperatures and efficiency; while response time is enhanced by increasing the air mass flow 
rate (decrease in NJ, one would not want to maximize air mass flow rate in most designs 
because by doing so the air experiences a negligible rise in temperature. 

The effect of the length of the heat exchanger is examined with the N w parameter. Figs. 9 
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and 10 show the outlet water and air temperatures, respectively versus the dimensionless time 
with N w as a parameter. The x-axis is normalized with respect to the flush time of the water 
of the N w = 1 case so that time durations are equivalent between the different N w curves. Thus, 
one flush time is shown at Tau’=l for the N w = 1 case, while for the N w =8 case, the flush time 
occurs at Tau’=8 since its length is 8 times longer than for N w =l (assuming the same fluid 
velocity in each case). As N w decreases, the outlet water temperature approaches the inlet water 
temperature, the outlet air temperature approaches the water temperature, and steady state is 
reached more quickly than for the larger values of N w . The transient efficiency curves shown 
in Fig. 1 1 more directly compare response time for the N w variations. Note that N w = 1 ,4, and 
8 correspond to the conventional units of NTU=0.5, 2.0, and 4.0, respectively, here. 

Approximate solutions are available for the cases examined here. The current work is not 
intended to replace previously available solutions, but rather to complement them. Effort has 
been aimed at outlining and describing the results and limitations of a simplified numerical 
solution approach to solve for heat exchanger transients. 

Extension of the current work is ongoing to further utilize the flexibility of the numerical 
approach. Areas being studied include examining the effects of tube rows on transient behavior 
and examining changes in secondary fluid temperatures and the flow rates of both fluids. 
Comparisons of solutions between the present modeling approach and a simplified model using 
ordinary differential equations is also underway. It is anticipated that these efforts will add 
further insight into modeling requirements as well as insight into design of heat exchangers with 
consideration of transients. 
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NOMENCLATURE 


A Area (m 2 ) 

c p Specific heat (J/kg-K) 

C r Dimensionless ratio of metal capacitance to water capacitance. 

Q Capacitance rate (W/K), mass flow specific heat product, 
h Convective heat transfer coefficient (W/m 2 -K) 

Dimensionless ratio of water side resistance to air side resistance. 

L Total flow length experienced by water (m) 

LMTD Logarithmic mean temperature difference: (T M -TJ/ln[(T m -T u )/(T m -T M )] 
m Mass flow rate (kg/s) 

NTU Number of transfer units, UA/C^. 

N w (NTU») Dimensionless number, inverse of water side resistance-capacitance product. 

N, (NTUJ Dimensionless number, inverse of air side resistance-capacitance product. 

P Tube inner perimeter (m) 
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t Time (s) 

UA Overall heat transfer coefficient (W/K) 

V Velocity (m/s) 

d Volume (m 3 ) 

x Axial position along tube (m) 

greek 

p Density (kg/m 3 ). 

if Surface efficiency 

r Dimensionless time, IVJL. 

9 Dimensionless temperature (T-TJ/(T wi -TJ 

X Dimensionless position, x/L. 

subscripts 
w water 

a air 

ai inlet air 

ao outlet air 

i grid point 

m metal (fin plus tubing) 

f fin 

t tube 
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■with NTUw-1.0, NTUa- 1000.0. CR-1.0. Hr- 1.0 



Tau (time/flush lime) 


Fig. 1 . Comparison of numerical and exact solutions (dimensionless water outlet 
temperature versus dimensionless time) for the case of negligible air side 
capacitance, with the grid/time spacing ratio as a parameter. 
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with NTUw*1.0, NTUa- 1000.0, CR-1.0, Hr* 1.0 
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Fig. 2. Comparison of numerical and exact solutions (integrated dimensionless air outlet 
temperature versus dimensionless time) for the case of negligible air side capacitance, 
with the grid/time spacing ratio as a parameter. 
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Fig. 3. Dependency of outlet dimensionless water temperature versus dimensionless 
time on the metal/ water capacitance ratio, with all other parameters held fixed. 
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Fig. 4. Dependency of the dimensionless integrated air outlet temperature versus 
dimensionless time on the metal/ water capacitance ratio, with all other 

parameters held fixed. 


with NTUw-1.0, Hr- 1.0, NTUa-1.0 



Fig. 5. Dependency of the transient efficiency versus dimensionless time on the 
metal/water capacitance ratio, with all other parameters held fixed. 
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with NTUw-1.0, Hr- 1.0, Cr- 1.0 



Teu (time/flush lirne) 


Fig. 6. Dependency of the dimensionless water outlet temperature versus dimensionless 
time on the N, parameter, with all other parameters held fixed 
(simulated air capacitance rate effects). 


with NTUw-1.0. Hr- 1.0, Cr-1.0 



Fig. 7. Dependency of the integrated dimensionless air outlet temperature versus 
dimensionless time on the N, parameter, with all other parameters 
held fixed (simulated air capacitance rate effects). 
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■with NTUw-l.O, Hr- 1.0, Cr- 1.0 



Tau (time/flush time) 


Fig. 8. Dependency of the transient efficiency versus dimensionless time 
on the N, parameter, with all other parameters held fixed 
(simulated air capacitance rate effects). 
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Fig. 9. Dependency of the dimensionless water outiet temperature versus normalized 
dimensionless time on the N w parameter, with all other parameters 
held fixed (simulated heat exchanger length effects). 
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■ 

Tau (tim a /normalized flush time) 


Fig. 10. Dependency of the integrated dimensionless air outlet temperature versus 
normalized dimensionless time on the N w parameter, with all other parameters 
held fixed (simulated heat exchanger length effects). 



Fig. 11. Dependency of the transient efficiency versus normalized dimensionless time 
on the N w parameter, all other parameters held fixed 
(simulated heat exchanger length effects). 
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